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摘 要 : 以 室外 试验 为 基础 ,应 用 饱和 - 非 饱和 土壤 水 分 运 移 理论 ,利用 Hydrus 软件 对 甘肃 黑 方 台 包 气 带 水 分 运 
移 规律 进行 数值 模拟 ,分 析 土 壤 剖 面 不 同 深度 处 含水 量变 化 规律 ,将 模拟 结果 与 实测 结果 进行 对 比分 析 , 并 将 拟 合 
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到 两 种 不 同 的 介质 中 。 结 果 表 明 :Q 土壤 含水 率 呈 峰 和 谷 波动 变化 ;@) 均 质 黄土 层 和 上 覆 5 cm 厚 的 粉 质 


壤土 的 黄土 层 中 ,土壤 水 分 入 渗 率 和 潜水 藻 发 率 均 随 潜水 埋 深 的 增加 而 减少 ,在 潜水 埋 深 300 cm 以 下 时 ,两 种 介 
质 的 土壤 水 分 人 渗 率 趋 于 同一 值 , 含 粉 质 培土 的 潜水 蒸发 率 小 于 均 质 黄土 的 潜水 蒸发 率 。 表 明 该 含水 量 的 变化 具 
有 时 空 特征 ,土壤 水 分 人 渗 率 潜水 蒸发 率 随 着 潜水 埋 深 变化 ,其 规律 整体 相似 ,这 对 研究 滑坡 水 分 分 布 具有 重要 
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甘肃 省 永 靖 县 黑 方 台 自 20 世纪 60 年 代 以 来 ， 
先后 发 生 滑坡 116 次 ,累计 37 人 死亡 ,100 多 人 受 
伤 , 直 接 经 济 损失 达 1.7 x10 元 ,实际 生产 生活 中 ， 
长 期 大 水 漫灌 是 引起 黑 方 台 边 坡 失 稳 ,以 致 滑坡 反 
复发 生 的 最 主要 的 诱发 因素 "“”。 在 半 干 旱地 区 ， 
地 下 水 的 主要 补给 来 源 有 降水 灌溉 水 和 地 表 水 的 
AŻ, 地 下 水 的 消耗 途径 有 蒸发 等 ”5 。 尽 管 黄 土 具 
有 很 强 的 渗透 能 力 ,可 是 在 强烈 的 大 气 营 发 作用 下 ， 
有 限 的 降水 还 是 很 难 通过 巨 厚 的 包 气 带 层 而 大 面积 
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移 特征 与 植被 的 根系 密度 等 特征 有 关 , 并 且 与 土壤 
的 结构 及 其 构造 相关 性 比较 大 。 

本 文 利用 德国 公司 生产 的 Ku -pF 非 饱 和 导 水 
率 测定 系统 ,测定 相关 水 分 运 移 参数 ”" 。 针 对 植被 
稀少 的 荒地 的 水 分 运 移 规律 ,利用 Hydrus - 1D 软件 
构建 降雨 和 蒸发 条 件 下 黑 方 台 黄 土 层 包 气 带 水 分 运 
移 模 型 , 拟 合 土壤 水 分 运 移 参数 ,分 析 不 同 土壤 深度 
下 含水 量 的 变化 规律 ,并 对 比 上 和 覆 有 5 cm HURE 
黄土 的 土壤 水 分 入 渗 率 和 潜水 燕 发 率 随 潜水 埋 深 的 


入 渗 形 成 地 下 水 , 故 在 干旱 半 干 旱地 区 研究 包 气 
带 水 分 运 移 规律 至 关 重 要 。 研 究 土 壤 水 分 运 移 相关 
参数 的 方法 有 仪器 测量 等 的 直接 法 和 数值 模拟 等 的 
I8 BE 7r 35779 , 拟 合 土壤 水 分 特征 曲线 的 模型 有 
van Genuchten 方程 及 其 修正 方程 Brooks — Corey 方 
程 .Dual - porosity 77 f£ , Gardner Jy gi 等 , 其 中 
Gardner 方程 的 计算 简便 , van Cenuchten 方程 应 用 
广泛 ,再 通过 相关 软件 的 模拟 或 者 运用 MATLAB 编 
程 等 方法 来 获取 相关 参数 77 。 本 文选 取 van Ge- 
nuchten 方程 ,运用 Hydrus 软件 拟 合 土壤 水 分 特征 
曲线 的 参数 。 目 前 ,研究 干旱 地 区 有 植被 的 耕地 的 
土壤 水 分 运 移 多 采用 原 位 试验 、 同 位 素 或 者 数值 模 
拟 等 方法 "“-” 。 研 究 表明 ,有 植被 的 黄土 水 分 运 
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变化 特征 ""。 该 研究 成 果 对 黄土 水 分 运 移 规律 及 
黄土 滑坡 治理 具有 重要 价值 。 


1 材料 与 方法 


试验 区 概况 

黑 方 台 为 泽 水 河 与 黄河 交汇 处 的 四 级 阶地 台 
面 中 ,面积 13.5 km? 左右 , 台 缘 高 差 约 120 m™, 
多 年 平均 降水 量 187.6 ~317 mm , 丰 水 年 达 400 mm 
以 上 ,雨季 (7 一 9 月 ) 降 水 占 全 年 降水 量 的 70% ~ 
80% 。 试 验 区 在 黑 方 台 台 震 建 立 气象 站 ,收集 实时 
气象 数据 。 利 用 ECH,0 土壤 水 - 温度 监测 系统 , 监 
测 剖 面 温度 .含水 量 。 典 型 点 选择 在 斜坡 泉水 溢出 
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带 上 游 200 cm 部 位 ,垂直 方向 坦 人 ECHO 土壤 水 
温 监 测 系统 。 监 测 剂 面 深 200 cm ,土壤 剖面 每 隔 10 
em 布置 一 个 探头 ,如 图 1 所 示 。 


气温 监测 
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图 1 


典型 点 剖面 示意 图 


Fig.1 Schematic diagram of typical point profile. 


1.2 水 分 运 移 模型 的 建立 

模型 的 理论 基础 是 一 维 饱 和 非 饱 和 带 水 分 运 移 
Richards 77 f ^ ,和 Van Genuchten ?! 确定 的 土壤 
水 分 运动 曲线 公式 。 

1.2.1 水 分 运动 基本 方程 


| 


| 8-8, 
g 0, = 0, 


" (1) 


=[1+ (alhl)"]™” 
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K(8) =K0.[1-(1-0")"] 


式 中 :0 表示 含水 量 , em! «em 25: 表示 时 间 ,d;h 
表示 压力 水 头 ,cm;z 为 垂直 坐标 轴 , 将 坐标 原点 选 
在 地 面 , 向 上 设 为 正 ;0,,0.,0. 分 别 表示 饱和 含水 
量 .残余 含水 量 和 有 效 土 壤 含 水 量 ,， cm - 
n 表示 经 验 参 数 ( 曲 线性 状 参 数 ),a 单位 为 ( 
量 纲 ;!/ 表示 地 下 水 位 埋 深 ,cm;m =1 - 
1/n,0 <m <1;K(9) 表 示 渗 透 系 数 ,cm/d;K, 表示 饮 
和 渗透 系数 ,cm/d。 

1.2.2 边界 条 件 及 初始 条 件 ”在 原 位 试验 的 基础 
上 ,水 流 模型 的 上 边界 是 大 气 边界 ,为 第 二 类 边界 条 
件 , 输 入 实测 的 降雨 或 莹 发 数据 ;下 边界 在 潜水 面 
处 ,为 深层 排水 的 边界 条 件 。hCritA 为 地 表 人 允许 的 
最 大 负 压 值 ,采用 100 000 cm。 初 始 条 件 输入 野外 
条 件 下 实测 含水 量 。 


-3 
cm 


Q, 


cm™'),n 


0(z,t) =O (2), tot, 初始 条 件 


ðh 


-K(% -1 ) =q (t) ,上 边界 条 件 


q(h) = -Ayexp( B, Ih - GWLOLI) , 下 边界 条 件 

(2) 
式 中 :9 表示 初始 土壤 剖面 含水 量 , cm : ecm; 
qg (RR ERKE E ecm. h q(h) 表示 排水 速 
率 ,cm * h^ ;4,, .B,, 表 示 经 验 参 数 ,分 别 为 -0.1 和 
-0.02; GWLOL 表示 地 下 水 位 的 参考 位 置 ,200 cm, 
1.2.3 单元 划分 “研究 区 为 试验 介质 分 为 均 质 和 
EH 5 cm 厚 的 粉 质 壤 土 的 黄土 层 。 黄 土 层 含水 量 
变化 主要 发 生 在 0 ~200 cm 土壤 深度 之 间 529 ,模型 
将 0 ~200 cm 的 土 层 划 分 为 5 层 ,观测 点 设置 在 30、 
50,100,150 cm 深度 处 。 

空间 及 时 间 单 位 分 别 采用 cm、h。 初 始 时 间 步 

长 1b, 计算 过 程 中 根据 收敛 的 情况 ,时 间 步 长 由 程 
序 自动 调整 ,最 小 时 间 步 长 为 1 b, 最 大 时 间 步 长 为 
5h, 
1.2.4 土壤 参数 识别 及 模型 求解 ”基于 实测 的 土 
坏 含 水 率 与 室内 测定 的 土壤 水 分 特征 参数 ,利用 试 
估 - 校正 法 对 土壤 参数 进行 识别 汪 ) 。 土 壤 含水 率 
的 变化 对 a 不 敏感 ,a 取 值 在 0. 001 ~ 0. 005 之 间 ， 
而 nn 对 含水 率 的 变化 最 敏感 ;残余 含水 率 与 饱和 全 
水 率 的 变化 不 大 。 因 此 ,在 参数 识别 中 主要 对 n 进 
行 识别 。 

本 文 Hydrus - 1D 模型 只 考虑 水 分 运动 ,模型 中 
方程 求解 采用 Calerkin ” 线性 有 限 元 法 ,模拟 裸 地 
只 有 降水 和 蒸发 下 水 分 的 迁移 与 转化 过 程 。 

现今 人 们 已 经 发 展 了 室内 外 试验 手段 来 测定 水 
分 特征 曲线 ,但 大 多 数 方法 费时 费力 ,试验 结果 存在 
明显 的 不 确定 性 ,而 且 土 壤 在 水 平和 垂 向 的 空间 变 
异性 也 限制 了 直接 测定 方法 在 实际 工作 中 的 应 
用 5 。 本 次 试验 选取 原状 样 (1 号 和 2 号 ) 与 扰动 
样 (3 号 和 4 号 ) ,采用 ku - pF 非 饱和 导 水 率 测定 系 
统 ,自动 测定 非 饱 和 导 水 率 和 水 分 特征 曲线 ,应 用 
Matlab 最 小 二 乘法 对 试验 数据 进行 参数 拟 合 “。 
参数 见 表 1。 

利用 相关 系数 和 协 方差 对 拟 合 值 和 实测 值 之 间 
进行 误差 分 析 , 由 相关 系数 和 协 方差 可 以 看 出 拟 合 
值 和 实测 值 拟 合 效果 好 ,具体 见 表 2。 
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表 1 土壤 水 分 特征 参数 拟 合 表 


Tab.1 Fitted values of the soil water characteristic 


parameters 
ü 特征 参数 
号 介质 k, [A 0, e n 
/( em : d!) /(cm? em) /(em em 3) [om 
1 原状 样 14.184 0.109 0.365 — 3.98x10-? 4.54 
2 原状 样 19.728 0.099 0.3444 3.77 x10? 4.20 
3 扰动 样 26.064 0.091 0.306 | 2.97 x10 ? 2.58 
4 ji 22.022 0. 094 0.373  4.24x10? 3.82 
表 2 误差 分 析 表 
Tab.2 Error analysis 
编号 ”介质 ER scil E 
相关 系数 (7) 协 方差 (Cov) 
1 原状 样 0. 994 5.312 x10-3 
2 原状 样 0. 994 4.893 x 10 ? 
3 扰动 样 0. 989 4.748 x 10 2 
4 扰动 样 0. 994 5.880 x10 ? 


3 结果 


3.1 蒸发 降雨 下 土壤 含水 量变 化 过 程 

由 于 含水 量 分 布 受 时 空 变化 影响 ,将 研究 时 
段 划分 为 3 个 时 段 即 冻 期 (2012 — 12 - 10 一 2013 - 
02 - 10) 、 融 期 (2013 -02 - 10—2013 -05 - 10) ,dE 
冻 期 (2013 - 05 - 10—2013 - 10 - 10) ,根据 试验 数 
据 , 分 析 冻 期 和 融 期 每 隔 20 d , 非 冻 期 每 隔 40 d 的 
自 地 表 至 土壤 深度 200 cm 的 含水 量 的 变化 特征 。 

由 图 1 看 出 ,在 冻 期 ,黄土 层 含水 量 随 测 深 的 增 
大 而 增加 ,降雨 旦 入 渗 状 态 。 在 剖面 上 ,不 同时 期 土 
壤 深 度 30 cm 以 上 的 含水 量变 化 较 大 ,30 em 之 下 
相差 小 ,但 整体 变化 规律 相同 ,黄土 层 含水 量 随 土壤 
深度 的 增 大 而 增加 ,并 且 随 着 土壤 深度 的 增加 含水 
量 趋 于 相同 。 

由 图 2 看 出 ,在 融 期 ,相同 土壤 深度 随 着 时 间 的 
延长 含水 量 增加 。 土 壤 深 度 10 cm 附近 含水 量 相 
近 , 随 着 土壤 深度 增加 ,含水 量 整 体 呈 增 大 趋势 , 当 
深度 至 200 cm 时 黄土 层 中 含水 量 相同 。 

由 图 3 看 出 , 非 冻 期 含水 量变 化 随 土壤 深度 的 
增 大 而 增加 ,但 变化 趋势 不 明显 。 在 同一 土壤 深度 
下 含水 量 递 增 。 土 壤 深 度 40 cm 以 上 含水 量 分 布 极 
不 均匀 , 极 差 较 大 ,分 析 原 因 是 该 时 期 蒸发 强烈 ,对 
地 表 附 近 含 水 量 影 响 较 大 。 相 同 土壤 深度 条 件 下 含 
水 量 呈 先 增 大 后 减少 的 规律 ,因为 该 时 期 虽然 蒸发 
强烈 ,但 是 相 比 前 期 ,该 期 降雨 集中 ,所 以 ,强烈 蒸发 


的 表层 土壤 合 水 量 低 ,在 7 月 强 降 雨 时 ,表土 含水 量 
增加 ,接近 饱和 。 由 图 1 ~3 可 以 看 出 ,含水 率 呈 现 
峰 谷 波动 变化 ,不 同时 间 峰 谷 位 置 相同 。 在 30 em, 
50 cm 130 cm 170 cm 处 为 波峰 ,40 cm, 70 cm 140 
cm 处 为 波 谷 。 
3.2 ”模型 的 识别 与 验证 

将 2013 年 5 月 10 日 至 7 月 10 日 和 7 月 10 日 
至 10 月 8 日 设 为 识别 斯 和 验证 期 。 基 于 野外 实测 

含水 量 /(cm? * cm?) 

0.10 0.15 0.20 0.25 0.30 
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图 2 冻 期 黄土 剖面 含水 量 随时 间 变 化 曲线 


Fig.2 Changing curves of soil water content distribution 


with time in freezing weather 
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图 3 融 期 黄土 剖面 含水 量 随时 间 变 化 曲线 


Fig.3 Changing curves of soil water content 
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distribution with time in melting period 


ae -8 1 N = 
含水 量 (cm cm”) SD = JFE equ (3) 
0.15 0.25 0.35 0.45 
A i RP: Y, 表示 样本 模拟 值 ; Y, 表示 样本 实测 值 ; N 
20 
30 上 表示 观测 样本 数 ;i 表示 样本 号 。 
40 - 、 Z 
50 | m 由 表 3 可 以 看 出 ,不 同 土壤 深度 含水 量 的 SD 
Bud ET 均 小 于 8% ,可 以 认为 模拟 得 到 的 水 分 运 移 参数 合 
E | A 2013-07-30 理 , 从 而 得 到 黄土 的 水 分 运 移 参数 ( 表 4) o 
gh 100 r —o— 2013-09-08 m z . 
X 110 上 RITE 表 3 不 同 土壤 深度 含水 量 的 模拟 值 
dd 与 实测 值 的 误差 
140 上 Tab.3 SD of Simulated and measured values 
150 f in different depth 
160 上 
170 r 土壤 深度 /cm SD/% 
180 r 
bul 30 6.7 
200 - 50 5.32 
图 4 ， 非 冻 期 剖面 含水 量 随时 间 变 化 曲线 e 
150 7.27 


Fig.4 Changing curves of soil water content distribution 


with time in freeze — free period 


表 4 土壤 水 分 运 移 参 数 拟 合 结果 


数据 ,经 Hydrus - 1D 数值 模拟 绘制 图 5 ,通过 识别 Tab.4 Fitting values of soil hydraulic parameters 
期 ,将 7 月 10 日 实测 值 与 模拟 值 对 比 ,经 过 验证 期 ， 参数 均 质 黄土 
验证 10 H 8 日 含水 量 的 拟 合 效果 , 校 验 建立 的 数值 6,/ (cm? * cm?) 0.098 23 
e 模型 的 合理 性 。 并 比较 同一 深度 不 同时 间 含 水 量 的 0/(cm * em?) 0.435 
N 拟 合 效果 ,如 图 6 所 示 。 o/ (cm!) 3.019 x10 7 
SY — vw " cm- d^! 
模拟 值 与 实测 值 的 拟 合 程度 采用 标准 差 SD om pois 
(standard deviation ) 定量 表示 , 即 n 1.466 5 
V SK E/(cn? * cem?) 含水 量 /(cma * cm?) 
:一 0.00 0.10 0.20 (0.30 0.40 0.50 0.00 0.10 0.20 0.30 0.40 0.50 
C 10 -一 10 
20 上 。 20 上 
30 | (a) 07-10 30 上 
40 | (10-08 
39 - 实测 什 50 上 
60 - 60l 
70 - 5 , 70 上 
80 - PUE B s0} 
Eg 90r p: 90r EE 
i 100r EK 100} . 实测 什 
g 110 | mE 一 一 模拟 值 
E 120} 120 上 。 
4 130 上 130 上 。 
140 - 140 上 . 
150 上 150 上 e 
160 上 160 上 . 
170 上 170 上 
180 - 180 上 . 
190 - 190 上 . 
200 - 200 - . 


图 5 2013 年 7 月 10 日 和 10 月 8 日 含水 量 模拟 值 与 实测 值 对 比 曲 线 


Fig.5 Comparison between simulated and measured values of soil water content on July 10 and October 8. 
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不 同 土壤 深度 下 含水 量 模拟 值 与 实测 值 对 比 曲线 


Fig.6 Comparison between simulated and measured values of soil water content in different depth 


0.50 f (c) 深度 100 cm 
E 
o 
E 
Š 
nim 
* 
4 
0.00 
07-05 07-25 08-14 09-03 09-23 10-13 
日 期 
图 6 
表 5 两 种 介质 的 土壤 水 分 运 移 参 数 
Tab.5 Parameters of soil hydraulic parameters 
of two soils 
参数 均 质 黄土 粉 砂 壤土 
0,/ (cm? - cem ^?) 0.098 23 0.067 

0,/ (cm * cm?) 0.435 0.45 

o/ (cm ^! 0.003 02 0.02 

k (cm * d^! ) 22.50 10.80 

n 1.466 5 1.41 


3.3 ”模型 应 用 与 分 析 

经 识别 和 验证 ,研究 区 介质 参数 误差 满足 精度 
要 求 ,可 将 此 参数 应 用 在 试验 场地 或 类 似 区 域 的 黄 
土 包 气 带 模型 应 用 上 。 在 2013 年 7 月 10 日 至 10 
月 8 日 薰 发 降雨 条 件 下 ,模型 应 用 分 别 模拟 两 种 不 
同 的 介质 ,一 种 模型 介质 是 均 质 黄土 , 另 一 种 介质 是 
上 层 是 5 cm 厚 的 粉 砂 壤土 的 黄土 ,介质 的 水 分 特性 
参数 如 表 5 ,用 以 分 析 两 种 介质 下 的 土壤 水 分 人 渗 
KAKE Az RB REKER RE SUBE 
3.4 土壤 水 分 入 渗 率 

由 图 7 可 知 , 均 质 黄土 的 土壤 水 分 人 渗 率 ,在 洪 
水 埋 深 为 10 ~ 80 cm 时 , 随 着 潜水 埋 深 的 增加 ,土壤 
水 分 人 渗 率 减 小 , 埋 深 40 ~50 cm 的 土壤 水 分 人 渗 


率 从 2.482 cm - d'KZ 1.618 cm - d”, 埋 深 80 
cm 以 下 的 土壤 水 分 入 渗 率 趋 于 稳定 值 1. 618 em - 
d”; 相 比 之 下 上 覆 5 cm 粉 质 壤土 的 介质 在 埋 深 10 
~40 cm 的 土壤 水 分 人 渗 率 从 1.660 em d 增加 
到 1.748 cm - d”, 埋 深 40 ~300 cm 的 土壤 水 分 人 
渗 率 有 小 幅 的 波动 , 埋 深 300 cm 以 下 的 土壤 水 分 人 
渗 率 与 均 质 黄土 趋 于 同一 稳定 值 1.618 em - d, 
3.5 潜水 蒸发 率 

由 图 8 可 知 ,在 潜水 埋 深 为 10 ~70 cm , 均 质 黄 
土 的 蒸发 强度 随 着 潜水 埋 深 的 增加 而 减 小 , 埋 深 40 
~50 cm Z RIRE M 2.489 em - d! eI 38 1.652 


3.00 f 


一 一 均 质 。 一 * 一 含 夹层 
2.50 V. 

2.00F | 
1.50 F 


1.00 r 


土壤 入 渗 率 /(cm * d) 


0.50 r 
0.00 
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埋 深 /cm 
图 7 两 种 介质 结构 土壤 水 分 入 渗 率 模拟 值 


Fig.7 Fitting values of infiltration rate of two soils 
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图 8 两 种 介质 结构 土壤 潜水 蒸发 率 模 拟 值 


Fig.8 Fitting values of phreatic water evaporation 


rate of two soils 


cm * d^ , 埋 深 80 em D F BUZA Az REGES T fae fH 
1.68 em * d ; 相 比 之 下 ,上 覆 5 em 粉 质 壤土 的 介 
质 在 埋 深 10 ~40 cm 的 蒸发 强度 从 1.641 em. d^" 
增加 到 1.716 em . d^! , 埋 深 40 ~ 150 cm 的 蒸发 强 
度 有 小 幅 的 波动 , 埋 深 150 ~ 300 cm 的 蒸发 强度 从 
1.645 cm - d 趋 于 稳定 值 1.110 em - d ;整体 上 ， 
均 质 黄土 的 蒸发 强度 大 于 含 粉 质 壤 土 介 质 的 蒸发 强 
度 。 


5 讨论 


以 往 对 黄土 层 含水 量变 化 研究 深度 都 大 于 10 
m, 时 间 短 ,其 变化 规律 呈现 峰 谷 波动 “” 77s o 
文 研究 深度 为 200 em ,在 冻 期 , 随 着 土壤 深度 的 增 
加 含水 量 基 本 趋 于 一 致 ; 融 期 ,同一 土壤 深度 下 含水 


量 随时 间 延 长 呈 递 增 的 趋势 ,在 地 表 附 近 含 水 量 相 
近 , 当 土壤 深度 至 200 cm 时 黄土 层 中 含水 量 趋 于 一 
致 ; 非 冻 期 ,同一 土壤 深度 下 含水 量 递增 ,地 表 附 近 
因为 蒸发 强烈 ,地 温 变 化 较 大 以 至 于 售 水 量 分 布 极 
不 均匀 , 极 差 较 大 ,但 相同 土壤 深度 条 件 下 含水 量 呈 
先 增 大 后 减少 的 变化 趋势 ,而 随 着 土壤 深度 的 增加 
含水 量 趋 于 相同 ,这 可 能 与 黄土 的 土壤 质地 和 颗粒 
组 成 有 关系 …”-”。 本 文 可 以 看 去 ,在 冻 期 含水 量 
分 布 比较 稳定 , 随 着 时 间 的 延长 ,气候 的 变化 ,气温 
上 升 ,表层 土壤 解冻 ,水 分 下 渗 , 随 着 深度 的 增加 , 含 
水 率 呈 现 峰 谷 波动 ,不同 时间 峰 谷 位 置 相同 。 在 30 
cm,50 cm 130 cm 170 cm 处 为 波峰 ,40 cm, 70 em, 
140 cm 处 为 波 谷 。 

目前 大 多 数 研 究 土 壤 质 地 、 含 水 量 以 及 时 间 与 
潜水 蔡 发 率 的 关系 “两 种 介质 结构 的 土壤 水 
分 人 渗 率 和 洪水 燕 发 率 与 潜水 埋 深 有 关 。 上 和 覆 有 5 


cm 厚 的 粉 质 壤土 的 介质 的 土壤 水 分 入 渗 率 随 潜 水 
埋 深 变化 的 趋势 比 均 质 黄土 的 平缓 ; 埋 深 300 cm 以 
下 时 ,两 种 介质 的 土壤 水 分 人 渗 率 趋 于 同一 稳定 值 
1.618 cm * d [BÆ EE 5 cm 厚 的 粉 质 壤土 介 
质 的 潜水 莹 发 率 小 于 均 质 黄土 的 潜水 蒸发 率 , 与 冻 
融 期 的 变化 规律 相似 ””) 。 整 体 趋势 都 是 先 降低 最 
后 趋 于 稳定 值 ,可 以 看 出 ,在 细 颗 粒 上 禾 粗 颗粒 的 黄 
土 介 质 时 ,可 以 降低 土壤 水 分 人 渗 率 和 潜水 蒸发 率 ， 
有 效 地 保护 水 分 的 流失 ,但 同时 也 阻止 了 水 分 的 人 
渗 补 给 。 

在 治理 滑坡 时 ,往往 忽略 水 分 运 移 特 征 , 实 际 
上 ,水 分 运 移 分 布 对 滑坡 的 形成 及 产生 具有 重要 的 
作用 。 模 拟 并 分 析 水 分 在 黄土 中 的 运 移 规律 ,对 滑 
坡 的 预防 具有 一 定 的 积极 作用 。 模 型 的 缺点 是 未 考 
虑 水 分 的 二 维 运动 ,试验 结果 存在 一 定 的 误差 ;而 且 
缺少 具有 植被 土壤 的 水 分 运 移 特 征 ,研究 介质 比较 
单一 。 本 文 研究 黄土 包 气 带 水 分 的 一 维 运动 规律 ， 
该 方法 针 省 时 省 力 具 有 很 好 的 拟 合 效果 。 针 对 室内 
小 型 试验 ,也 可 以 进一步 设计 方案 来 验证 拟 合 效 果 
的 稳定 性 ,所 以 ,可 以 把 HYDRUS -1D 应 用 在 分 析 
黄土 土壤 水 分 运 移 规 律 的 研究 上 。 

以 往 对 土壤 水 分 人 渗 率 和 潜水 蒸发 率 的 计算 多 
采用 经 验 公 式 估算 ”3 , 上 且 研 究 对 象 多 以 农田 林地 
的 植物 生长 或 灌溉 为 主 , 但 是 黄土 地 区 地 质地 貌 类 
型 复杂 ,面积 广阔 ,黄土 是 岩 土 颗粒 水 .空气 三 者 同 
时 存在 的 一 个 复杂 系统 ,人 研究 少 且 理论 基础 不 完善 ， 
而 采用 经 验 公式 进行 估算 缺乏 代表 性 , 存在 较 大 误 
差 “"。 本 文 在 研究 区 选取 典型 点 设计 原 位 试验 ， 
通过 实测 数据 的 模拟 ,分 析 黄 土 包 气 带 水 分 运 移 规 
律 。 
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«d - 
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在 研究 区 原 位 试验 的 基础 上 ,对 黄土 包 气 带 降 
水 入 渗 率 与 土壤 蒸发 强度 进行 数值 模拟 ,由 观测 数 
据 和 计算 结果 分 析 得 到 以 下 结论 : 

(1) 黄土 层 包 气 带 含 水 量变 化 具有 时 空 特点 。 
空间 上 随 着 深度 的 增加 , 呈 峰 谷 波动 变化 ;时 间 上 ， 
在 同一 深度 随 着 时 间 的 延长 含水 量 先 增 加 后 减 小 ， 
且 到 200 cm 时 含水 量 趋 于 一 致 。 

(2) 模拟 计算 了 两 种 岩 性 组 合 的 黄土 包 气 带 降 
水 入 渗 率 与 土壤 蒸发 强度 ,得 到 黄土 包 气 带 降 水 人 
渗 率 和 土壤 莱 发 强度 受 包 气 带 岩 性 结构 和 水 位 埋 深 
的 影响 。(3 ) 拟 合 得 到 的 黄土 水 分 运 移 参数 ,可 以 
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应 用 到 该 区 的 水 文 地 质 研 究 中 。 
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Simulation of loess soil water migration in certain northwest 


Abstract: According to analyze the saturated — unsaturated soil hydraulic parameters, vadose water migration 
mechanism is numerically simulated in the research by using Hydrus - 1 D software , which was based on the outdoor 
test and saturated. — unsaturated soil moisture migration theoretical basis . The model has been verified using the 
field experiments ,and the fitting parameters are applied to two soils, homogeneous soil lager and overlying 5 cm 
thick silty loam soil lager. The results showed that 1) the moisture content was fluctuating in peak valley. 2) Soil 
water infiltration rate and the phreatic water evaporation rate decreases with an increasing of the phreatic water level 
in two soils; two kinds of soil water infiltration rate tend to be the same value when the phreatic water level was less 
than 300 cm , medium containing silty loam of the evaporation rate is less than the phreatic water evaporation rate of 
homogeneous loess. In conclusion,the change of water content has spatial and temporal characteristics. The rate of 
two soils water infiltration and the phreatic water evaporation are similar with the change of phreatic water level. It 
provides scientific basis for the water distribution and change of the landslide. 


Key words:  Loess; The vadose zone; Water migration; Numerical Simulation; Water migration mechanism 


